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Abstract 

We have proposed to use an effective theory to describe interactions of an TV TV-system. The 
effective theory can be constructed in analogy to the existing effective theory for an TViV-system. 
In this work we study the next-to-leading order correction to NN scattering near the threshold 
in the effective theory. We find that the experimental data can be well described with the 
effective theory. 

Interactions between a nucleon and an antinucleon have been studied extensively with potential 
models(see [H [2j |3] and references therein). In general the potential of these models includes 
various effects of meson exchanges and is supplemented with a phenomenologically determined 
imaginary part to account for the NN annihilation. The descriptions of the low energy NN 
scattering with these potential models are successful. Recently there is a renewed interest to study 
NN interactions, stimulated by experimental observations of the threshold enhancement of a pp 
system in the radiative decay of J/^>[3], -B-decays[5j and e + e _ -annihilations[6]. It has been shown 
that the enhancement can be explained with the final state interaction, where the enhancement is 
determined by the NN scattering amplitude near the threshold [71 [HI l9l fTTl [T2l [13l fl^l fT5l fT6l fT7] . 

Because the momentum transfer in NN scattering near the threshold is small, the inner struc- 
ture of a nucleon or antinucleon can not be seen. This indicates that one can take the nucleon and 
antinucleon as point-like particles and construct an effective theory to describe the NN scattering 
near the threshold. Near the threshold the interactions through 7r-exchanges are most important 
among those of exchanges of other mesons, hence, the effective theory only contains N-, N- and 
7r-fields. In [17] we have proposed such an effective theory. We have used the tree-level result of 
the effective theory and a partial sum of higher order effects to successfully describe the observed 
enhancement and the NN scattering near the threshold. In this work we study the next-to-leading 
order correction in the effective theory. 

An effective theory for an NN system near the threshold can be constructed in analogy to the 
effective theory for an NN system. The effective theory for the iViV system has been proposed 
in[l8l [T9| [20] and studied extensively in [I9j [20l [21]. In constructing such an effective theory one 
makes an power expansion in the momentum near the threshold. In the effective theory, the U.V. 
divergences are regularized with the dimensional regularization. The interactions with 7r-mesons 
are fixed with chiral symmetry. A power counting to determine the relative importance of different 
terms in the effective theory has to be established. But, there are distinct difference between the 
effective theory of NN interactions and that of NN interactions. 

It is well-known that for an NN system the scattering lengthes are large and there is a shallow 
bound state in 3 Si — 3 Z?i channel and a virtual state in 1 Sq channel. To take these facts into account 
in the effective theory, the power subtraction scheme has been introduced [20J. With this scheme a 



systematic power counting of the effective theory can be established. In the case of NN systems, 
the scattering lengthes, according LEAR experiment[7] and model results [E], are around lfm. They 
are much smaller than those of NN systems. Therefor one can use the minimal subtraction scheme 
and hence the simple power counting for the effective theory of NN systems. 

Another difference is that an NN system can be annihilated into mesons, while a NN system 
can not be annihilated. The annihilation of an NN system into virtual or real mesons results in 
that the dispersive and absorptive part of the NN scattering amplitudes are relatively of the same 
importance. In order to incorporate this fact some coupling constants in the effective theory of an 
NN system are complex numbers. One should keep in mind that the complex coupling constants 
here do not mean the violation of time-reversal symmetry. The effective Lagrangian with complex 
coupling constants should be understood as for the purpose to effectively build the S*-operator 
and hence scattering amplitudes. This can be understood as the following: One can imagine that 
the effective theory is obtained from a perturbative matching of a more fundamental theory. In 
the more fundamental theory with the time-reversal symmetry scattering amplitudes can become 
complex beyond tree-level because absorptive parts can exist at one- or more loop level. The 
imaginary parts of the coupling constants in the effective theory are from these absorptive parts in 
the matching. If the underlying theory respects to the time-reversal symmetry and other discrete 
symmetry like charge-conjugation and parity, the generated operators through the matching in the 
effective theory can not violate these discrete symmetries of the underlying theory. This fact also 
tell us that the effective theory should be built with these operators which are C-, P- and T-even. 

At low energy, i.e., for an NN system near its threshold, we need to consider theory contains 
contact interactions and interactions with pions which is consistent with chiral symmetry. At low 
energy we can use the nonrelativistic fields to describe the nucleon N. The nonrelativistic fields 
are given as 

*-(£)■ *=(£)■ « 

The field ip{x ) annihilates a nucleon(an antinucleon) and the field ip'(x) creates a nucleon(an 
antinucleon) . The pion- field is defined with Pauli matrices t 1 acting in the isospin space as: 

with / ff « 93 MeV. The field ^{i = 1,2,3) is real. Under SU(2) L x SU(2) R chiral symmetry the 
fields transform as 

£ -> L£U ] = U£R\ £ = f -> LER\ V -»• Uij), (3) 

where L, R are global transformations in SU(2)l and SU(2)r respectively and U is a pion-field- 
dependent SU(2) matrix. From £ we can give out the vector and axial-vector pion currents as 

V, = 1(^+^0,0, A^-W^-tfdrf), (4) 

where the axial current and the chiral covariant derivative = (d^ + V^) transform linearly 
as 

UAprt, Dp -> UD^UK (5) 
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With the above fields the effective Lagrangian can be constructed as: 



C eff 
C-N 



^(iD + D 2 /2m)^ + xHiDo - D 2 /2m) X + ji^hxH + jV^WtV 

+^VxxW + ^VV'xxVrty + W f 0- ' + g A X ] <r ■ Ax + C7(p 2 ) 

f 2 m 2 f 2 

^Tr^E^X + -^Tr(£ + £t) + C(p 4 ), 



r 



(6) 



where 5,4 ~ 1.25. cj with the isospin index I is the coupling constant in the Sq channel, while dj 
is the coupling constant in the 3 Si channel. These coupling constants are in general complex to 
account for the NN annihilation. 

The power counting of the effective theory is discussed in detail in[17j. It is easy from the 
effective Lagrangian to derive the tree-level NN scattering amplitude. The amplitude is an expan- 
sion in the three momentum p of N or N. The leading order is at order of p°. the amplitude at 
this order is determined by those contact terms given explicitly in Eq.(6) and the interaction from 
exchange of single tt. We will consider the momentum region p ~ m n . From Eq.(6), the NNir- or 
NNtt vertex is proportional to p. In the single 7r-exchange, the power of p from the two vertices is 
canceled by the power of the denominator of the 7r-propagator. This is why the amplitude with sin- 
gle 7r-exchange is taken at order of jr. The next-to-leading order is at order of p and is determined 
by one-loop diagrams. It should be noted that some of one-loop diagrams will not contribute at 
the order. They can contribute at order of p 2 or higher. As we will show that it is easy to find 
those one-loop diagrams which give the contributions at order of p. At order of p 2 the amplitude 
receives corrections from different sources. The corrections can come from those one-loop diagrams 
which contribute at order of p 2 and come from two-loop diagrams. They also come from those 
contact terms in the effective Lagrangian at order of p 2 , which we have not included in Eq.(6). 
Those contact terms in general contain derivatives. We will use the dimensional regularization and 
the minimal subtraction scheme as discussed before. 



Figure 1: The leading order diagrams. Fig. la is for the contact terms, Fig. lb is from one ir- 
exchange. 

After giving the effective theory, we now consider the NN scattering near the threshold: 



of nucleons. As discussed before, the leading order contribution is at 0(jP). The contributions at 
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this order can be represented by diagrams in Fig.l. The leading order amplitude can be expressed 
as, 



7"(o) 



mt 



(8) 



where £ and rj are the spinors of the nucleon and antinucleon respectively. 

At one-loop level, not all diagrams give contributions at the order 0(p 1 ). At one- loop level one 
meets a loop integral of a four momentum k. One can perform the fc°-integration first by using a 
contour in the complex A;°-plan. The contribution of the integration comes from poles insider the 
contour. The poles can appear from denominators of nucleon propagators or 7r-propagators. It is 
easy to show that only those contributions from the poles of nucleon propagators are at the order 
0(p 1 ), while other contributions are at higher orders. With this in mind, only those diagrams at 
one-loop given in Fig. 2 give contributions at the next-to-leading order. 





(b) (c) — tar 

Figure 2: The diagrams which give the next-to-leading order corrections. 



It is rather straightforward to evaluate these loop contributions. In fact results of some loop 
integrals exist in the literature. Hence, we give here the results without showing here detailed 
calculations. It should be noted that all contributions in Fig. 2 are finite in the dimensional regular- 
ization. As discussed in the above, we will use the minimal subtraction scheme for renormalization. 
In this case, the next-to-leading order contributions do not depend on the renormalization scale 
fi, or the subtraction point. The /U-dependence of coupling constants starts at order of p 2 . The 
contributions from Fig. 2a read: 



T 



2<i 



.mk 

' 4-7T 



2 2 j2 j2 



(9) 



r 

T 



yr( c )£ . I (k)^a i T a T^a i T a r/ + h(k)^a ■ kr a T^a ■ kr a 



we note that this contribution with the minimal subtraction scheme can be equivalently calculated 
by taking a cut cutting the nucleon loop in Fig. 2a. 
The contributions from Fig.2b and Fig. 2c read: 

■ 1 / 9 A ■ 
1 t 9A 

2c 

In the above r( c ) is a matrix acting in the spin- and isospin space. The sum of the direct product 
r( c ) <g> r( c ) is implied. The sum reads: 



-i—{—) 2 ^r^7] ■ I {k)rfa i T a T^a i T a £ + h(k)rfa ■ pr a T^a ■ pr a £ 
47r 2f n L 



(10) 



r (c) gr (c) = £0j^ T , C l,a^,a , ^0 _j „ , dl 



4 4 



a H — -a t Qs> cr t . 
4 



(11) 
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The functions Iq and 1\ are from loop-integrals. They are: 
h(k) 



mk 



l + y-iy^2y~ + y(y + 2) (~ ln(l + ^) + itan 1 y ^ 
1 - 3y + 3% - ?/ 2 + 3y) — ln(l + - + itan -1 J- 

V 2 y V y 

2 



2/c 2 



(12) 



The contributions from Fig. 2d are the most difficult parts. However, the relevant results of loop 
integrals can be found in [22J. We have numerically checked these results and found an agreement. 
The contributions can be written in the form: 



r 



2,1 



F v 5 %3 + F w (k - p) l (k - P y + F T (k x p)\k x p)i 
+iF P (k x pf feV^r/ - ^irfoS + 2£Vr a £? ? t T a r ? - 2^T a ^T a a i 7 ] ) } , (13) 



the functions Fsy,T,W,P are given as: 

1 



F S 
F w 
F P 
F T 



{ ^ {4(2m 2 + q 2 )T (p) + 2q 2 T l (p) - 4ip - (2m 2 + q 2 ) 2 G (p, ?)} 

— \f v = ^-G 2 (p,q), 
q z Att 

i- (2r (p) + 2T l (p) - (2ml + I 2 ) [Go(p, q) + 2G?i(p, g)]} , 
1 

47T 



{G (p,?)+4£i(p,?) + 4G 3 (p, g )}, 



where the complex- valued functions are given as, 

r (p) 



2p it 
arctan hdn — 



yVi 2 + 4p 2 , 



ri(p) 

G (p, q) 

Gi(p,q) 
G 2 (p,q) 
Gs(p,q) 



,., [m n + ip- (ml + 2p 2 )T (p)] 



■ Qm-n . pq + R 
arcsm \- i in 



urn.,, 

,2 , 



i? = m4 + p 2 w 2 , w = \J 4m 2 + g 2 , 



1 

2^ 
1 
2p 
1 

r (p)-2^( g )-(m 2 +2 P 2 )G (p, g ) 1 , 

V^? ' A[q) ~ Yq arCtan 2^ 

/G (p,(z) + (mj + 2p 2 )G 1 (p,q) +A(q), 
lTi(p)-p 2 G (p,q) - 2(m\ + 2p 2 )G 1 (p,q) 

4^2 _ ^2 



where we define q = \k — p\ here. The total results can be written as: 

T(p, k) = T (0) (p, k) + T m (p, k) + 0(p 2 ) 1 



(14) 



(15) 



(16) 
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where is given in Eq.([8]) and the order 0(p 1 ) result is: 



T (i) = T 


+ r 


+ T 


+ T 






2a 


2b 


2c 


2d 



(17) 



With the above results we have the complete result for NN scattering near the threshold at the 
next-to-leading order. We will use this result to fit the experimental data. 




Figure 3: The narrow and long bubble represents the amplitude Aya'gj] (E) at certain order with 
the subtraction of contributions from physical cuts or cut diagrams. The sum of diagrams can be 
performed. 



Our results for the scattering amplitude can be further improved by including some higher 
order corrections without a concrete calculation. We have noticed in [161 [T7] that some corrections 
from higher orders can be summed into a compact form. The summation can be done for partial 
waves. We define our partial waves from the NN scattering amplitude Ti(p, k, s±, s 2 , s^, s 2 ) with 
the isospin I = 0, 1 as: 

Ti(p,k,s 1 ,s 2 ,s' 1 ,s 2 ) = ^Yee 3 (jp/p)Y t , e t 3 (k/k)(l/2, si; 1/2, s 2 |s, s 3 )(l/2, si; 1/2, s' 2 \s' ', s' 3 



■?,) 



<i, is, s, s 3 \j,j 3 )(i', 4 s', s' 3 \j\j' 3 )5 jr 5 jd ,J ss ,T bu/StI] (E). (18) 

In the above the repeated indexes are summed and E is the total kinetic energy E = p 2 /m of 
the system. The summation can be explained with Fig. 3. We will take s-waves as examples to 
illustrate this. 

Supposing we have calculated the scattering amplitude 7joooo,o] at certain orders and we 
denote this amplitude as -4.[oooo,o](-^)) i n which possible contributions from physical cuts or cut 
diagrams are subtracted. At tree-level, -Apooo.o] (E) is just the tree-level amplitude. At the next- 
to-leading order -4[oooo,o] (E) is determined by diagrams from Fig.l and those from Fig.2 after 
subtracting the contributions of cut diagrams of Fig.2. In Fig. 3 we denote this amplitude with 
the narrow long bubble. Now using this amplitude .Aroooo.o] (E) one can generate those diagrams 
as ladder diagrams at higher orders as shown in Fig. 3. Physically the interpretation of Fig. 3 is 
the following: The iViV undergoes a multiple scattering process AW — > NN —)••••—> AW. Each 
scattering is due to the amplitude «4.[oooo,o] (-^0 Each pair of NN is on-shell in Fig.3. We will call 
amplitudes for such a multiple scattering process as rescattering amplitudes. It is easy to shown 
that the sum of these diagrams is the sum of a geometric series. One then finds the sum for the 
amplitude as: 



7[oooo,o] {E) — ^[oooo,o] (E) 



1 ~ i J^2 ^[oooo,o] (E) 



+ •■•• (19) 



This can be generalized to other partial waves. We will improve our result by adding the rescattering 
effects. The improvement is done only for S- waves with the assumption that the amplitudes with 
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small I are dominant because of the finite interaction range. A small complication is with the case 
of s = 1 is that there is a mixing of /S-waves and D-waves. In this case, the summation takes a 
form of 2 x 2 matrix. Details can be found in |17| . 



t> 
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0.05 
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Figure 4: The fitting results for the cross-section of pp elastic scattering. The dashed line is the 
leading order fit. The solid line is the next-to- leading order fit. 

We have performed numerical fits with our results for experimental data. The data are taken 
from [23J. In |17j we have studied the pp scattering and the enhancement in J/ip — > "fpp and 
e + e~ — > pp together with a combined fit. For the enhancement, certain assumptions have been 
made in [T7] to relate the observed enhancement to NN scattering near the threshold. In fact, the 
study of the enhancement can also be done with the approach of effective theories in a consistent 
way combined with the effective theory here for NN system. We leave this to future works. Here 
we only focus on pp elastic scattering in order to see how good our effective theory works. In Fig. 4. 
we give our fitting result with the amplitude at leading and next-to-leading order, respectively. 
With our results we are able to describe the experimental data below E ~ 50MeV, corresponding 
to p ~ 200MeV. The fitting quality with the leading and the next-to-leading result is given by 
x/d.o.f. w 2.6 and x/d-o.f. sa 1.7, respectively. At tree-level we can not determine the coupling 
co,i or do t i separately, because they appear as the combination as cn + ci or do+di in the amplitude. 
With the NLO results we can determine these coupling constants in the unit of GeV -2 as: 



where the fitting errors are given in the brackets. For c\ the errors are very small in comparison with 
the numbers. We hence do not give the errors. From Fig. 4. one can see that our NLO amplitude 



(180±4)+i(20±3), ci = -211-M87, 

-(60 ± 2) + i(61 ±4), di = (37 ± 3) + i(56 ± 4). 



(20) 



7 



can describe experimental data well. With the determined couplings we can also determine the 
scattering lengthes of S- waves: 

a [1So 0] = -(5.88 ± 0.22) - i(1.12 ± 0.17), ap S(J)1 ] = 0.67 - i0.09, 

ap Slj0] = 4.24 - i0.06, ap 5l)1] = -(0.25 ± 0.06) - i(1.08 ± 0.08), (21) 

in unit of fm. In the above the last integer in [• • • ] indicates the isospin / = or I = 1. 

We have also performed the numerical fits by adding the rescattering effects into the leading- 
and next-to-leading amplitude as discussed in the above. The fitting results are represented in 
Fig. 5. From Fig. 5. one can realize that with the rescattering effects the experimental data can also 
be well described, indicated by x/d-o.f. pa 1.7. From the fit with next-to-leading amplitude and 
rescattering effects we have for the coupling constants in unit of GeV -2 as: 

c = (70 ± 37) + i(0 ± 190), ci = (132 ± 92) + i(0 ± 248), 

do = -(25±2)+i(0±14), d x = -(38 ± 8) + i(99 ± 80). (22) 

The corresponding scattering lengthes are 

op Soj0] = (0.3 ± 2) - i(0 ± 10), a[i So ,i] = (0.34 ± 0.09) - i(0.00 ± 0.24), 

a[ 3 Si 0] = 4.21 - »(0 ± 0.01), op Sl) i] = (1.20 ± 0.15) - i(1.91 ± 1.54). (23) 

The determined scattering lengthes are qualitatively in agreement with results from some models [8j 
and from LEAR experiment [7]. They are smaller or much smaller than those of an AW system. 
This fact supports the argument for using the minimal subtraction scheme and hence the standard 
power counting for our effective theory. 

Although the data can be well described from Fig. 5, the significant changes in some coupling 
constants indicate that corrections from higher orders can be significant. In Fig. 5 we also give 
results for the contributions from those partial waves with I ^ where the contribution from the 
partial wave 3 D X is subtracted. These contributions essentially come from 7r-exchanges. From Fig. 5. 
we can see that these contributions will become dominant when the total energy E is approaching 
to the value around 50MeV, corresponding to p ~ 200MeV. This indicates that the expansion 
in p of the effective theory converges slowly or applicable region of the effective theory is small. 
This has also been found in the effective theory of an AW system. The study in |21j shows that 
the contributions from 7r-exchanges become dominant and much larger than predictions at leading 
orders around p ~ lOOMeV in some spin triplet channels. The reason for this is the true expansion 
parameter for 7r-exchanges is at order of (g\m w m) / '{Stt f%) . Because the nucleon mass m is large 
in comparison with f n , the expansion parameter is not small. In our case in which contributions 
from 7r-exchanges become dominant at p ~ lOOMeV, the situation seems better than that for 
AW systems, because the coupling constants of the local interactions, which give the important 
contributions at leading order, are roughly twice larger than the corresponding coupling constants 
in the effective theory of an AW-system. 

It is interesting to compare the approach of the effective theory with the approach of the the 
effective range expansion for describing AW scattering. In the effective range expansion, e.g., in 
|24j . the scattering amplitude, or more exactly, the phase-shift is expanded in p and only first two 
terms in the expansion are kept. The first term is the scattering length and the second term is 
determined by the interaction range. These terms can be complex. If we neglect interactions with 
7r in our effective theory, our effective theory will only consist of local operators. The perturbative 
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Figure 5: The fitting result for the cross-section of pp elastic scattering. The dashed line is the 
leading order fit. The solid line is the next to leading order fit while including the rescattering effect 
in 1 So and 3 Si — S D\ channel. The dotted line is the contribution from 1 So and 3 Si — 3 Di with 
the rescattering amplitude. The dashed-dotted line is the contribution from all the other partial 
waves up to the next-to-leading order. 



expansion of coupling constants of the local operators can be summed into a compact form. In 
the case, our approach will be identical to that of the effective range expansion. However, the 
existence of interactions with tt makes two approaches different. We notice that interactions with 
7r can not be represented with local operators. To use local operators for the interactions, one in 
essence makes an expansion in p/m n for the contribution from Fig. lb. But this expansion converges 
very slowly with the convergence range about several MeV, reflecting the fact that the scattering 
amplitude through 7r-exchanges has a cut very near- and below the threshold. In our approach we 
do not make such an expansion in p/m n . 

To summarize: We have studied an effective theory for an NN system near the threshold. 
The effective theory is constructed by an expansion of momenta near the threshold and including 
interactions with 7r-mesons fixed with chiral symmetry. The power counting of the effective theory 
has been given. We have obtained the next-to-leading order results for NN scattering and compared 
with experimental data. The data can be well fitted with our predictions. But, there can be the 
problem that the expansion of our effective theory does not converge quickly. This may need to be 
studied further. If there are more experimental data near the threshold, our effective theory can 
be tested more accurately. We notice that one can also predict the spin-dependent NN scattering 
by using our effective theory and compare with results from the planned experiment [25]. This will 
provide another interesting test of our effective theory. 
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